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INTRODUCTION 


The  first  signal,  i.e.,  the  P  wave,  arriving  at.  a  seismograph  from  a  shallow,  tolese- 
ismio  event  will  he  mixed  with  its  own  reflection  off  the  surface  of  the  earth.  The 
recorded  signal,  s(t  ) ,  may  be  represented  as  the  P  wave,  p(  t).  of  the  event,  plus  its 
undisl.ort.ed  echo  reversed  in  phase  and  delayed  by  time  d: 

s{t)  -  p(t)  -  r  •pit  -d)  (1) 

where  r  is  the  reflection  coefficient,  which  can  range  from  0  to  1.  If  !\u)  is  the  power 

spectrum  of  p,  then  the  spectrum  of  the  signal  s  is  (Cohen,  1970): 

P(u)  *  [  1  +  r2  -  2*r*cos  (u*  ri,)] 

If  the  reflectivity,  r,  is  1,  as  it  tends  to  be  at.  the  lower  frequencies,  then  the  spectrum 
of  the  recorded  signal  is  the  spectrum  of  the  P  wave  shaped  by  the  factor 

2\  1  -  cos  (u*  d)\ 

Thus  the  signal  spectrum  will  go  to  zero  at  every  frequency,  a  ,  at.  which 

u  =  2*n*n/d  (n  =  i  0, 1,<J,D.„.)  (^) 

If  the  power  spectrum  of  the  P  wave  itself  is  flat,  then  the  spectrum  of  the  recorded 

signal  will  have  a  sinusoidal  shape.  There  will  be  a  series  of  regularly  spaced  holes  cen¬ 
tered  on  frequencies  the  locations  of  which  will  depend  on  d  ,  the  delay.  That  is,  the 
hole  locations  will  depend  on  d  expect  for  the  hole  centered  on  frequency  zero,  which 
is,  of  course,  always  at  the  one  location. 

The  factor  shaping  the  signal  amplitude  spectrum,  i.e.,  the  square  root  of  the 
power  spectrum,  is 

V2  *|  l-cos(  u  *  d)  ]°  5 

'Hie  size  of  this  hole  produced  in  the  amplitude  spectrum  by  the  interference  of  the 
reflected  signal  may  be  characterized  by  its  half  width,  i.e.,  its  width  where  it  has 
reduced  the  spectrum  by  one  half.  This  occurs  when 


~  V  V.V'.^  V 


<"  v-  «.* 


u  =  (2*  n*  Tikn /3) /d  (3) 

As  an  example,  the  size  of  the  holes  produced  in  an  amplitude  spectrum  by  an  echo 

delayed  0.55  seconds  is  0.6  hertz.  That  is.  they  will  have  a  half-width  of  0.6  hertz.  In 
general,  the  shallower  the  event  the  wider  the  spectral  hole. 

The  low-frequency  hole  in  the  P-wave  spectrum  of  shallow  seismic  events  is  a 
potentially  powerful  discriminant  for  distinguishing  the  seismic  signals  of  earthquakes 
from  those  of  underground  nuclear  explosions.  When  the  hole  is  diagnostic  of  a  depth 
below  the  reach  of  present  technology  the  event  must  a  natural  earthquake.  Some  of 
the  advantages  of  the  low-frequency  hole  as  a  discriminant  arc: 

a)  At  low  frequencies  the  coefficient  of  reflectivity  approaches  1.0,  ensuring  that  the 
spectral  hole  is  not  obscured  by  low  reflectivity  at  the  earth's  surface,  as  it  may  be  for 
the  fundamental  mode  and  it.s  harmonics,  at  some  test  sites. 


b)  Unlike  the  fundamental  mode  of  pP  interference  and  its  harmonics,  the  com¬ 
ponent  at  the  lower  extreme  has  an  unequivocal  location  in  the  spectrum.  There  can 
be  confusion  about  the  position,  and  therefore  about  the  identity,  of  the  other  holes 
produced  by  interference.  The  hole  of  the  fundamental  mode  in  the  L0NGSH0T  spec¬ 
trum.  for  instance,  was  anomalously  located  (Cohen.  1970).  Hut  the  very  low  frequen¬ 
cy  component,  because  of  its  unique  position  cannot  be  confused  with  spurious  holes 
in  the  spectrum,  such  as  those  produced  by  discrete  multipaths.  (An  example  of  such 
multipatfiing  is  found  in  the  record  of  0NCSH0T  at  Kskdalemuir.  Scotland,  which  re¬ 
veals  two  discrete  direct  ray-paths  between  the  source  and  the  receiver;  Douglas  et 
al..  19  t'd.)  I  hough  multipathing  can  cause  spectral  scalloping,  it  does  not  produce  a 
hole  at  the  lower  extreme  of  the  spectrum. 
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c)  Scattering  of  surface-waves  to  body-waves  in  the  vicinity  of  the  event  is  reduced  at. 
the  lower  frequencies  since  the  scatterers  in  the  crust  are  smaller  than  l  he 
corresponding  wave-lengths.  This  ensures  that  energy  which  has  not  undergone  the 
pP  cancellation  does  not  fill  in  the  spectral  hole  and  conceal  it,  as  it  may  do  at  the  fun¬ 
damental  frequency  and  its  harmonics. 

d)  It  cannot  be  disguised.  The  flat  spectrum  characteristic  of  deeper  events  cannot 
be  simulated  as  can  a  discriminant  such  as  complexity,  which  characterizes  earth¬ 
quakes  but  which,  conceptually  at  least,  can  be  imitated  by  an  array  of  nuclear  shots 
fired  in  rapid  sequence. 

The  difficulties  associated  with  the  use  of  the  very  low  frequency  hole  as  a  discrim¬ 
inant  are  these: 

1)  The  low  frequency  part  of  the  signal  is  not  optimumly  recorded.  The  short-period 

reeording  systems  in  the  field  have  responses  which  drop  off  sharply  at  the  low  end  of 
t  "’m.  The  long-period  systems  operate  just  in  the  frequency  range  of  interest 

but  at  law  gan. :  that  permit  the  observation  of  the  P  waves  of  only  larger  events. 

2)  The  six-second  microseisms  dominate  precisely  the  portion  of  the  spectrum  of  in¬ 
terest. 

3)  Though  there  are  now  mid-period  recording  systems  going  into  the  field,  they  are 
not  installed  in  arrays  and  thus  their  data  are  not  susceptible  of  the  analytic  tech¬ 
niques  required  to  overcome  the  microseisrn  noise,  which  is  well  organized. 


The  problems  associated  with  the  very  low  frequency  discriminant  are  not,  how¬ 
ever,  necessarily  insuperable.  Since  the  six-second  microseisms  are  well  organized 
they  are  susceptible  to  array  processing,  and  may  thus  be  removable.  There  are  abun¬ 
dant  data,  both  past  and  present,  from  short-period  arrays.  These  records  will  admit 
of  array  processing  which  may  allow  both  the  suppression  of  the  microseisms  and  the 
reconstitution  of  the  low-frequency  end  of  the  signal  spectrum. 


OBSKIfVATlONS  OF  Tl  IK  pi*  VKIfY  1.0W  r’SQUKNCY  IIOI.K 
IN  SlIOltT-l’KUIOD  KXIM.OSION  Sl'KCTKA 


In  the  present,  research  t.o  evaluate  the  feasibility  of  using  the  very  low  frequency 
hole  as  a  discriminant  we  first  examine  the  spectra  of  events  large  enough  to 
overwhelm  the  microseism  noise.  We  do  this  simply  to  verify  that  there  is  such  a  hole 
to  be  observed  in  the  short-period  data,  given  the  steep  roll-ofT  at  the  low-frequency 
end  of  the  response  of  the  short-period  systems.  We  choose  the  1.ASA  short-period 
subarray  records  of  the  Amchilka  test  1.0NGSI10T.  The  event  is  suitable  because  l)  it  is 
at  teleseismie  distance,  2)  it  is  large  enough  to  be  just  within  the  dynamic  range  of  the 
system,  thus  affording  the  maximum  signal-Lo-noise  ratio,  and,  U)  at  the  time  of 
LONGSHOT  the  entire  array  was  still  in  place,  having  2b  short-period  vertical  seismome¬ 
ters  per  subarray.  Also,  a  control  event  is  available  for  comparison.  An  earthquake 
took  place  a  few  months  afterward,  on  December  12,  19(>f>,  within  about  a  degree  of 
IjONGSHOT  and  O.bb  magnitude  unit  smaller,  and  the  1,ASA  recording  of  it  has  survived. 
The  estimated  depth  of  the  earthquake  is  49  kilometers. 

An  idea  oT  the  appearance  of  the  spectral  hole  we  may  expect  to  sec  is  given  in 
Figure  1,  which  shows  the  anticipated  effect,  on  an  otherwise  flat  amplitude  spectrum, 
of  pP  interference  for  the  depth  of  burial  of  LONGSHOT.  The  expression  for  this  spec¬ 
trum  was  derived  in  the  introduction.  A  coefficient  of  reflectivity  of  1.0  is  assumed;  the 
event  is  at  teleseismie  distance.  The  amplitude  spectrum  of  I.ONGSf iOT  should  have  a 
hole  of  this  size,  shape  and  location  superposed  on  it. 

However,  the  signal  as  recorded  has  been  reshaped  by  earth  response  and  instru¬ 
ment  response,  which  will  distort  the  features  we  are  attempting  to  identify.  Neverthe¬ 
less,  the  control  earthquake  has  also  been  through  the  same  physical  filters  and  its 
spectrum  can  be  used  as  a  template  against  which  to  compare  the  shot  spectrum  to 
locate  and  identify  its  distinctive  features.  Since  the  anticipated  source  spectrum  of 
an  earthquake  is  flat  at  the  lower  end,  up  to  the  corner  frequency,  we  can  expect 
1.0NGSH0T  to  have  an  amplitude  spectrum  which  differs  in  shape  from  the  spectrum  of 
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Figure  1.  The  effect  on  a  flat  signal  spectrum  of  a  reflection 
delayed  0.55  sec  after  the  signal  itself,  with  reflectivity  equal 
to  l.O.  (0.55  sec  is  the  observed  pP  delay  for  LONGSHOT;  Cohen, 
1970).  The  actual  LONGSHOT  spectrum,  shown  in  Figures  2-9,  rolls 
off  toward  zero  frequency  at  about  the  same  rate,  with  respect  to 
the  control  earthquake. 
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Uh:  control  earthquake  as  Figure  1  differs  from  constant,  amplitude. 

For  that  comparison  Figures  2  through  9  juxtapose  actual  I ’-wave  spectra  of 
LONGSHOT  with  those  of  the  control  earthquake.  F.ach  figure  contains  the  l’-wave 
bearnsum  spectrum  of  the  shot  superposed  on  that  of  the  earthquake,  both  from 
records  at  a  specified  UtSA  subarray.  The  two  events  were  of  different  magnitude  and 
so  to  facilitate  the  comparison  of  the  spectral  shapes  the  LONGSHOT  spectra  in  these 
figures  have  been  adjusted  vertically  such  that  the  maximum  amplitude  of  the  shot 
coincides  with  that  of  the  earthquake.  Thus  the  scale  on  the  ordinate  applies  only  to 
the  earthquake.  (The  adjustment  is  entirely  graphical  to  facilitate  visual  comparison; 
the  shot  data  themselves  were  not  normalized  to  the  earthquakes’.) 

We  use  subarray  beamsurns  rather  than  single  traces  for  computing  thesi^-SpCctra. 
In  so  doing,  we  maximize  the  signal-to-noise  ratios,  but  we  also  avail  ourselves  of  sta¬ 
tistical  control  over  the  signal  estimates  by  means  of  F  spectra,  which  are  defined  for 
bearnsum  spectra  but  not  for  those  of  single  traces,  as  discussed  below.  The  usual  limi¬ 
tation  on  beamforming  that  a  trace  can  only  be  shifted  an  integral  number  of  the  data 
sample-intervals  is  avoided  here.  To  form  the  beam  we  shift  a  given  trace  the  integral 
number  of  sample  intervals  nearest  to  the  exact  time-delay  required,  and  then  we 
Fourier  transform  that  record  and  shift  the  remaining  fraction  of  a  sample  interval  in 
the  frequency  domain.  Thus  we  avoid  possible  artifacts  in  the  spectra  that  might  result 
from  beaming  inaccuracies,  especially  at  higher  frequencies. 

In  each  of  these  figures  the  spectrum  of  the  LONGSHOT  P  wave  can  be  seen  to  fall 
below  the  earthquake  spectrum  uniformly  and  steadily,  from  the  peak  at  about  1.0  Hz 
back  toward  zero  frequency.  Between  that  peak  and  0.15  Hz  the  shot  spectrum  drops 
ofl  about  -14db,  on  the  average,  with  respect  to  the  earthquake.  Note  that  this  is  equal 
to  the  anticipated  fall  in  amplitude,  shown  in  Figure  1,  where  from  the  peak  amplitude 
at  1.0  Hz  the  spectrum  drops  to  about  -13db  at  0.15  Hz.  This  is  the  observation  we 
sought  to  make.  The  very  low  frequency  discriminant  survived  the  filtering  effects  of 
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re  3.  From  LASA  subarray  B3  beamsums:  the  P  spectrum  of  LONGSHOT  superposed 
that  of  the  nearby  control  earthquake  of  December  12,  1965.  Note  the  roll- 
of  LONGSHOT  toward  zero  frequency,  with  respect  to  the  quake,  due  to  pP  in- 
erence.  Refer  to  Figure  2. 
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off  of  LONGSHOT  toward  zero  frequency,  with  respect  to  the  quake,  due  to  pP  in 
terference.  Refer  to  Figure  2. 
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e  9.  From  LASA  subarray  F3  beamsums:  the  P  spectrum  of  LONGSHOT  superpose 
hat  of  the  nearby  control  earthquake  of  December  12,  1965.  Note  the  roll 
f  LONGSHOT  toward  zero  frequency,  with  respect  to  the  quake,  due  to  pP  in 
rence.  Refer  to  Figure  2. 


transmission  through  the  earth  and  through  the  sensor,  in  spite  of  being  only  margi¬ 
nally  within  the  range  of  the  seismometer  frequency  response.  Note  that  in  each  of 
these  records  in  Figures  2  through  9,  both  for  the  shot  and  for  the  earthquake,  the 
noise  is  at  least  -<?Odb  below  the  signal,  as  shown  in  the  example  in  Figure  10. 
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Figure  10.  The  LASA  AO-subarray  beamaum  spectra  for  the  LONGSHOT 
P  wave  and  the  preceding  noise.  The  plus  symbols  represent  the 
P-wave  spectrum. 


AN  ANOMALOUS  PEAK  AT  ZERO  FREQUENCY  IN  LASA  SHORT-PERIOD  SPECTRA 


We  would  like  to  be  able  to  identify  and  measure  the  very  low  frequency  discrim¬ 
inant  in  the  spectra  of  underground  nuclear  explosions  without  having  to  resort  to  the 
spectrum  of  a  control  earthquake  as  a  template  with  which  to  compare  them.  That 
necessity  imposes  an  extra  constraint  on  the  use  of  the  discriminant.  One  objection  is 
that  a  suitable  nearby  earthquake  may  riot  be  available.  It  is  desireable  to  be  able  to 
remove  the  instrument  response  and  estimated  earth  response  from  the  shot  record  to 
observe  and  measure  the  discriminant  hole  directly.  However,  in  the  case  of  the  LASA- 
LONGSHOT  data,  upon  correcting  for  the  amplitude  response  of  sensor  and  recorder  we 
find  instead  of  a  hole,  a  spike,  or  hump,  more  than  40db  above  the  rest  of  the  spec¬ 
trum.  The  significance  and  the  source  of  the  spike  are  not  immediately  apparent.  It  is 
present  not  only  in  the  shot  spectrum,  as.  for  example,  in  Figure  1 1;  it  is  present  in  the 
earthquake  spectrum,  Figure  12;  in  the  spectrum  of  the  beam  residual,  Figure  I'd;  and 
in  the  noise  spectrum,  Figure  14.  In  all  the  spectra  the  spike  is  approximately  propor¬ 
tional  to  the  rest,  of  the  spectrum,  i.e.,  if  is  present  in  the  spectrum  before  the  arrival 
of  the  signal,  and  when  the  signal  arrives  it  increases  proportionally. 

Though  it  is  not  at  once  clear  how  to  account  for  the  observed  hump  around  zero 
frequency,  it  is  possible  to  decide  whether  it  is  part  of  the  signal.  We  can  do  that 
because  our  spectra  are  computed  from  beamsurnmed  seismograms.  An  F  spectrum  is 
defined  for  every  beamsurri,  and  whether  a  band  of  energy  in  a  certain  frequency  range 
belongs  to  the  beam  may  be  inferred  from  that  F  spectrum.  The  F  spectrum  is  simply 
the  F  statistic  as  a  function  of  frequency,  arid  the  1’’  statistic  at  a  given  frequency  is  the 
ratio  of  the  energy  in  the  beam  to  the  residual  energy  at  that  frequency,  all  times  the 
number  of  channels  minus  1.  For  a  discussion  of  the  1’  statistic  in  signal  processing, 
see  Shorn  way,  1971. 

Figure  lb  shows  an  example  of  one  of  the  I.ASA-I.ONGSI !0T  short-period  beamsum 
spectra  and  its  accompanying  F  spectrum.  Note  that,  the  holes  due  to  pP  interference 
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Figure  11.  The  LASA  D4-subarray  beamsum  spectrum  of  L0NGSH0T 
corrected  for  instrument  response.  Observe  the  large  peak  at 
zero  frequency.  Similar  anomalous  peaks  appear  in  the  spectra  of 
the  control  earthquake,  the  shot  beamsum  residuals,  and  the  pre¬ 
shot  noise,  scaled  in  each  case  to  the  rest  of  the  spectrum,  as 
shown  in  Figures  12,  13,  and  14.  The  window  length  is  51.2  sec. 
There  is  no  smoothing. 
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Figure  12.  The  LASA  Cl-subarray  beamsum  spectrum  of  the  control 
earthquake  corrected  for  instrument  response.  Observe  the  peak 
at  zero  frequency.  Similar  anomalous  peaks  appear  in  the  spectra 
of  LONGSHOT,  the  shot  beamsum  residuals,  and  the  pre-shot  noise, 
scaled  in  each  case  to  the  rest  of  the  spectrum,  as  shown  in  Fig¬ 
ures  11,  13,  and  14.  The  window  length  is  51.2  sec.  There  is  no 
smo  o  t  h ing  . 
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Figure  1A.  The  LASA  AO-subarray  beamsum  spectrum  of  pr e-LONGSHOT 
noise  corrected  for  instrument  response.  Observe  the  peak  at 
zero  frequency.  Similar  anomalous  peaks  appear  in  the  spectra  of 
LONGSHOT  itself,  the  control  earthquake,  and  the  shot  beamsum 
residuals,  scaled  in  each  case  to  the  rest  of  the  spectrum,  as 
shown  in  Figures  11,  12,  and  13.  The  window  length  iB  51.2  sec. 
There  is  no  smoothing. 


at  very  low  frequency  and  at  the  fundamental  frequency  are  both  represented  in  the  F 
spectrum.  The  hole  at  the  fundamental  frequency,  around  1.9  Hz,  can  be  accounted  for 
by  the  drop  to  near  zero  of  the  amplitude  level  at  that  point  in  the  signal  spectrum, 
since  the  F  spectrum  is  proportional  to  the  estimated  signal-to-noise  ratio.  But  there 
is  a  hole  in  the  F  spectrum  at  very  low  frequency  as  well,  in  spite  of  the  amplitude  max¬ 
imum  at  that  point.  This  is  accounted  for  by  the  residual  energy,  that  is,  the  energy 
not  in  the  beam,  i.e.,  noise,  which,  notwithstanding  the  magnitude  of  the  energy  in  the 
signal  spectrum  at  this  point,  is  so  much  the  greater  and  reduces  the  signal-to-noise 
estimate  to  near  zero. 

The  significance  of  the  low  F  statistics  in  the  vicinity  of  zero  frequency  is  that  it  is 
improbable  that  much  of  the  energy  in  that  band  of  the  estimated  signal  spectrum  is, 
in  reality,  part  of  the  signal.  So  the  spectral  estimate  of  the  explosion  F-wave  indicates 
large  amplitude  near  zero  frequency,  but  its  F  spectrum  shows  that  that  energy  prob¬ 
ably  doesn’t  belong  to  the  signal  after  all.  Notwithstanding,  it  has  obscured  the  very 
low  frequency  spectral  hole  we  seek. 

Thus  the  large  amplitudes  near  zero  frequency,  which  are  orders  of  magnitude 
above  the  rest  of  the  beamsum  spectrum,  do  not  belong  to  signal,  that  is,  they  arc  not 
organized  and  are  not  part  of  the  plane  wave,  and  yet  they  appear  when  the  signal 
appears.  Nevertheless,  the  same  phenomenon  is  present  in  the  earthquake  records, 
and  so  in  lieu  of  explicitly  removing  this  and  other  extraneous  effects,  i.e.,  instrument 
and  earth  responses,  as  we  should  prefer  to  do.  for  the  present  we  can  continue  to  use 
control  earthquake  spectra  with  which  to  compare  the  explosion  spectra.  In  effect  we 
are  subtracting  the  log  spectrum  of  the  control  earthquake  from  that  of  the  shot., 
which  mounts  to  dividing  out  the  extraneous  effects  implicitly.  Still,  it  remains  a 
significant  objective  to  contrive  to  do  < il I  that  explicitly,  without  the  control  event, 
since  a  satisfactory  control  earthquake  may  not  exist 
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We  account  for  the  anomaly  of  large  amplitudes  near  zero  frequency  by  observing 
that  t  tic  bcamsum  is  a  biased  estimator,  notwithstanding  its  being  the  optimum  esti¬ 
mator  for  plane  waves.  Thus  beamsumming  reduces  (but  does  not  clinm;..te) 
incoherent  noise  by  a  factor  of  one  over  the  square  root  of  the  number  of  channels,  but 
the  noise  remnant  remains.  Hut  that  observation  suggests  that  the  remnant  itself  may 
be  estimated  and  removed  from  the  spectral  estimate,  correcting  the  bias,  and  indeed 
it  can.  In  the  Appendix  we  develop  a  function  to  apply  to  bcamsum  spectra  which 
corrects  for  the  bias.  The  development  is  based  upon,  and  the  bias  correction  factor  is 
a  function  of,  the  F  statistic.  The  factor  is  (F-l)/F,  and,  like  F.  it  varies  with  frequency. 

Applied  to  the  spectra  of  the  LONGSHOT  F  wave,  the  bias  correction  indeed 
removes  the  very  low  frequency  hump  in  the  l.ASA  bcamsum  spectra  of  subarrays  111 
and  F4,  as  shown  in  Figures  1(5  and  17.  (Note  that  the  only  other  correction  applied  to 
these  spectra  is  for  the  response  of  the  seismometer  and  recording  system.)  However, 
the  technique  doesn't  prove  effective  for  all  of  the  1ASA  subarray  beams.  Only  I  d  of 
them  respond.  In  the  case  of  the  others  the  F  statistic  around  zero  frequency  is  simi¬ 
larly  low,  but  unfortunately  not  sufficiently  low  to  make  the  correction  factor  small 
enough  to  effectively  suppress  the  dominant  amplitudes  in  that  band.  The  difficulty 
remains. 

Note  that  care  is  taken  in  the  spectral  computations  to  ensure  that  each  window  of 
data  has  zero  mean  and  zero  linear  trend,  and  a  window  contour  is  selected  to  minim¬ 
ize  side-lobes.  The  particular  window  chosen  is  one  of  several  optimum  windows  intro¬ 
duced  by  Nuttall  (1981). 

To  test  the  hypothesis  that  the  very  low  frequency  hump  is  an  artifact  of  the 
LASA-NORSAR  short-period  recording  system,  we  recovered  the  simultaneous  record¬ 
ings  of  LONGSHOT  made  on  the  long-period  seismometers  at  LASA,  which  were  collo¬ 
cated  with  the  short-period  subarray  center  elements.  The  hypothesis  can  be  at  once 
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substantiated  (or  definitively  eliminated)  by  comparison  of  the  spectra  of  simultane¬ 
ous,  identically  windowed  records  of  the  collocated  seismometer  pairs.  The  com¬ 
parison  may  also  ofTer  a  remedy  for  the  problem.  Since  the  response  of  the  short- 
period  system  at  low  frequency  is  not  based  on  actual  field  or  laboratory  measure¬ 
ments  but  is  merely  extrapolated  from  measurements  made  in  the  short-period  band, 
and  is  thus  open  to  question,  spectral  comparisons  offer  an  alternative  calibration 
technique  directly  employing  the  actual  response  measurements  made  in  that  fre¬ 
quency  band  on  the  long-period  system.  Should  such  calibration  of  the  very  low  end  of 
the  short-period  system’s  response  be  confirmed  by  repeated  measurements  and  com¬ 
parisons  with  the  long-period  data,  the  new  measurements  would  be  used  to  replace 
the  extrapolated  calibrations. 
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CONCLUSIONS  AND  UKCOM M KN OATIONS 


Wo  have  confirmat  ion  that  the  very  low  frequency  discriminant,  i.e.,  I  he  hole  .il 
very  low  frequency  due  to  pi1  interference  in  underground  nuclear  explosion  spec 
t.rii.  is  preserved  in  short-period  seismic-  recordings.  This  is  shown  rcpe.it.edly,  Iroie 
subarray  to  subarray,  in  comparisons  of  LONCSliOT  with  a  control  earthquake,  in 
recordings  of  those  events  made  at.  1.ASA.  This  means  that  a  powerful  discriminant, 
previously  neglected,  is  available  with  present  recording  systems.  It  remains  to  be 
seen  whether,  in  application,  the  discriminant  is  to  be  observed  directly  in  the  ex¬ 
plosion  spectra,  or  will  be  measured  in  comparison  with  spectra  of  control  earth¬ 
quakes,  as  in  this  report. 


Microseisms  occupy  and  obscure  precisely  the  spectral  band  in  which  we  wish 
to  observe  the  very  low  frequency  discriminant,  iiut  both  they  and  the  signals  of  in¬ 
terest  are  highly  organized  and  well  suited  to  maximum-likelihood  array  processing 
which  permits  the  unobstructed  observation  of  one  signal  in  the  presence  of  anoth¬ 
er.  We  have  developed  a  linear  high-resolution  frequency-wavenumber  process 
■which  operates  at  the  limit  of  such  refinements  and  is  eminently  suited  to  this  task, 
i.e.,  the  resolution  of  signals  not  well  separated  in  azimuth  and  velocity.  (Smart, 
1976).  In  its  linearity  it  has  an  advantage  not  shared  by  conventional  high- 
resolution  techniques:  its  spectral  estimates  are  undistorted  by  non-linear  opera¬ 
tions.  Spectral  distortion  is  a  significant  arid  undesirable  side  effect  of  non-linear 
manipulation. 

A  demonstration  of  this  capacity  to  recover  the  low-frequency  portion  of  the 
spectra  of  smaller  signals  in  the  presence  of  microseisms,  and  to  reveal  the  very 
low  frequency  discriminant  in  explosion  seismograms  is  urgently  called  for. 
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AITKNDIX 


The  correction  for  bcamsum  spectral  bias 


The  expected  value  of  amplitude  reduction 

1 


waveforms  from  an  array  of  seismometers  is 


VN 


achieved  by  summing  up  A’  nrr-y 
(It  is  assumed  that,  the  correlation 


between  waveforms  is  zero.  The  amplitudes  of  the  summed  waveform  are  normalized 
by  division  by  A'.) 

Thus,  notwithstanding  that  it  is  an  optimum  signal  estimator,  beam  summing  is 

biased.  It  doesn’t  altogether  remove  the  noise;  it  reduces  it  by  a  factor  of  — —■■f  (at. 

V  /V 

best).  The  remaining  noise  contaminates  resulting  spectral  estimates  of  signal.  How¬ 
ever,  that  residual  noise  (which  may  be  larger  than  the  signal  in  portions  of  the  spec¬ 
trum)  can  itself  be  estimated  from  the  I'-statistic.  Then  an  unuiasuu  signal  sp>  i  i  i.m 
may  be  computed. 


To  calculate  a  correction  factor  for  the  biased  spectrum  wo  compute  the  tin- 
expected  bias  (in  the  frequency  domain)  assuming  zero  correlation  between  the  noise 
at  one  station  and  that  at  another,  and  between  the  signal  and  the  the  noise  every¬ 
where.  The  beam-shifted  and  Fourier  transformed  array-recording  of  an  event  may  be 
represented  by  ,  where  d  is  the  complex  transform  coefficient,  at  a  given  frequency, 
of  a  single  channel  of  data,  and  n  is  the  channel  index  of  the  nF1  seismometer  station. 

Let  the  data  be  composed  of  S,  the  signal  itself,  plus  b„,  the  background  (noise). 

<**  =  S  +  bn  (A-l) 


Since  the  records  are  beam-shifted  the  delay,  or  step-out,  across  the  array  has  been 
removed  from  the  signal,  so  we  take  S  to  be  identical  from  channel  to  channel.  Then 
the  complex  Fourier  coefficient  of  the  (normalized)  beam  sum  at  the  given  frequency  is 
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77-2  ^  =  7 r£S+b» 


(Note  that  Fourier  transforming  and  summing  is  identical  to  summing  in  the  time 
domain  and  then  transforming.)  Let  the  subscripts  r  and  i  identify  the  real  and  ima¬ 
ginary  terms  in  the  complex  spectrum.  Then  the  squared  magnitude  of  the  coefficient 


jr£s+b» 


=  ^TNSr  +  £bnr  +  NSi  +  t  bni  2 
=  mz+  jf£bn  +  jj\sr£b„  +  st£bM 

71  {  n  n 


Note  that  the  cross  terms  in  'Tj~£  have  an  expectation 


of  zero  by  our  assumption 


that  the  noise  correlation  between  sensors  is  zero.  Then,  since  the  expected  value  of 
the  third  term  in  (3).  above,  is  also  zero,  by  our  assumption  of  zero  correlation 
between  signal  and  noise,  the  beam-sum  spectral  estimate  (squared)  has  the  expected 


value 


4-2  I0» 


at  each  frequency. 

Similarly,  the  expectation  of 


42  i  i'  =  42 1  ®  ♦  <>. 
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In  t  he  frequency  domain  t  he  l''-sl  al  istie  for  beam  sums  is  defined  as 


(A'-n  -Xx 

_ A  n  | 

.  W  i  ,  N 


Substituting  our  signal  model 


(N-l)  |S  |2+  — £  ,  i 2 
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Thus,  with  the  bias  removed,  our  estimate  of  the  signal  amplitude  squared,  at  a  given 


frequeney,  is 


Since 


Ibnl2 
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That  is,  the  beam  sum  amplitude  spectrum  times  the  square  root  of 
(F-l)/F  is  an  unbiased  signal  estimate.  {Note  that  ( F-\)/  F  is.  of  course,  fre¬ 
quency  dependent.) 

Note  that  the  range  of  expected  values  of  F  is  from  1  {  when 
|  5  |  goes  to  zero  )  to  ®  (  when  ^  \  bn  | 2  approaches  zero  ).  Notwithstanding  this 

n 

range  of  expected  values,  estimates  of  F  can  fall  below  one,  in  which  case  the  correc¬ 
tion  factor  becomes  ( F  —  1  )  /  F  <0.  Since  a  negative  amplitude  spectrum  is  mean¬ 
ingless,  we  take  the  factor  to  be  zero  when  F  <  1  . 
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